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Abstract. We study effects of turbulent mixing on the random growth of an interface 
in the problem of the deposition of a substance on a substrate. The growth is modelled 
by the well-known Kardar-Parisi-Zhang model. The turbulent advecting velocity field 
is modelled by the Kraichnan’s rapid-change ensemble: Gaussian statistics with the 
correlation function (vv) oc 6{t — t') where k is the wave number and 0 < ^ < 2 

is a free parameter. Effects of compressibility of the fluid are studied. Using the field 
theoretic renormalization group we show that, depending on the relation between the 
exponent ^ and the spatial dimension d, the system reveals different types of large- 
scale, long-time asymptotic behaviour, associated with four possible fixed points of the 
renormalization group equations. In addition to known regimes (ordinary diffusion, 
ordinary growth process, and passively advected scalar field), existence of a new 
nonequilibrium universality class is established. Practical calculations of the fixed 
point coordinates, their regions of stability and critical dimensions are calculated to 
the first order of the double expansion in ^ and e = 2 — d (one-loop approximation). It 
turns out that for incompressible fluid, the most realistic values ^ = 4/3 or 2 and d = 1 
or 2 correspond to the case of passive scalar field, when the nonlinearity of the KPZ 
model is irrelevant and the interface growth is completely determined by the turbulent 
transfer. If the compressibility becomes strong enough, the crossover in the critical 
behaviour occurs, and these values of d and ^ fall into the region of stability of the 
new regime, where the advection and the nonlinearity are both important. However, 
for this regime the coordinates of the fixed point lie in the unphysical region, so its 
physical interpretation remains an open problem. 


PACS numbers: 05.10.Cc, 05.70.Fh 


1. Introduction and description of the model 

Over decades, constant interest has been attracted to the growth processes in various 
physical systems: solidihcation and flame fronts, smoke and colloid aggregates, tumors, 
and so on; see e.g. n-i and references therein. A most prominent example is provided 
by the deposition of a substance on a substrate and growth of the corresponding phase 
boundary (interface). A number of microscopic models were proposed to describe those 
phenomena: Eden model |6], Edwards-Wilkinson model [7], restricted solid-on-solid 
model [8], ballistic deposition [9]; not an exhaustive list. 
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It turns out, however, that the growth processes share some important features with 
those of equilibrium nearly-critical systems: namely, self-similar (scaling) behaviour 
(power-like dependencies) with rather universal (independent of the details of a specihc 
process) exponents. In particular, the n-th order structure functions of a growth process 
behave as n-ia 

S„(i,r) = ([ft((,x) - ft(0,O)]"> ~ r = |x|. (1.1) 

Here h{x) = h{t,x.) is the height of the interface prohle, the brackets (...) denote 
averaging over the statistical ensemble, x cind ^ are referred to as the roughness 
exponent and the dynamical exponent, respectively, and Fn{-) is a certain universal 
scaling function. The asymptotic behaviour fll.ip takes place in the infrared (IR) range, 
where the time and space differences t, r are large in comparison with characteristic 
microscopic scales. 

It is then natural to try to describe universal properties of the growth processes on 
the base of a certain simplihed model for a smoothed (coarse-grained) height held, in 
analogy with the theory of critical state, where most typical universality classes (types 
of critical behaviour) are described by the classical <y9^-model [TITl ITT] . As the coarse¬ 
grained model of growth, one usually choses the Kardar-Parisi-Zhang (KPZ) model 
na, described by the nonlinear stochastic differential equation 

dth = M:od‘^h + Xo{dhy/2 + f. (1.2) 

Here the height held h{x) = h(t, x) depends on the d-dimensional substrate coordinate 
X, dt = d/dt, di = d/dxi, = didi is the Laplace operator and {dhY = dthdih; the 
summations over repeated tensor indices are always implied. The hrst term in the right- 
hand side of fll.2p describes the surface tension with the coefficient xq > 0. The second 
term represent an excess growth along the local normal to the surface. The parameter 
Ao can be of either sign; it can be scaled out, and in the following we set Aq = 1. 

Furthermore, / = f{x) is the Gaussian random noise with zero mean and given 
pair covariance 

{f{x)f{x')) = 2DoS{t - - x'), (1.3) 

with the positive amplitude factor Hq > oJil 

To be precise, the model (11.21) . (II.3p had appeared for the hrst time in the seminal 
paper by Forster, Nelson and Stephen [13] in terms of the purely longitudinal (solenoidal) 
vector held Ui = dih. Then, for Aq = —1 it represents the d-dimensional generalization 
of the Burgers equation. It can also be mapped onto a model of directed polymers in 
random media and on a model of Bose many-particle system with attraction; see e.g. 

[H]. 

Actually, the hrst two terms on the right-hand side of fll.2l) are just the simplest local 
ones that respect the symmetries h —)■ h-|-const and 0(d). Thus the KPZ model arises 

f Strictly speaking, a nonvanishing mean value (/) should be introduced in order to cancel a linear 
in-time growth of the mean value (h). Once we are interested in quantities like (11.11) that involve only 
differences of the fields, both the mean values can be simultaneously ignored. 
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ubiquitously in description of many nonequilibrium, disordered and driven diffusive 
systems. Then the held h can have different meanings. For example, in [15] the KPZ 
model and its ramihcations were used to study large-scale distribution of matter in the 
Universe. 

A few generalizations and modihcations of the original KPZ model were introduced: 
random noise with hnite correlation time [16], vector or matrix held Aim, modihed form 
of the nonlinearity [18] and anisotropic modihcations [19]. In connection to the latter, it 
is also worth to mention continuous anisotropic models of self-organized criticality IZD]. 

The powerful quantitative theory of the critical state is provided by the held 
theoretic renormalization group (RG); see the monographs [IDl El] and references 
therein. In the RG approach, possible universality classes are associated with IR 
attractive hxed points of renormalizable held theoretic models. 

The RG analysis of the KPZ model, pioneered in [T^ E2], eventually (after some 
misunderstanding) led to the following conclusions [21], [23]. The held theoretic version 
of the stochastic problem (ll.2p - (ll.3p is miltiplicatively renormalizable. The nonlinearity 
{dhY in fll.2p is IR irrelevant (in the sense of Wilson) for d > 2, logarithmic (marginal) 
for d = 2 and relevant for d <2. Thus it can be studied within the standard perturbative 
RG and the expansion in e = 2—d. The corresponding RG equations possess a nontrivial 
hxed point with the exponents y = 0, 2 ; = 2 (the exact relation y -|- .^ = 2 is dictated 
by Galilean symmetry). However, the hxed point for £ < 0 is IR repulsive, while for for 
£ > 0 it does not lie in the physical range of the model parameters {Dq, xq > 0) and 
thus can hardly describe the IR asymptotic behaviour of the problem. All these results 
are “perturbatively exact,” that is, exact in all orders of the expansion in e. 

One can nevertheless assume that the KPZ model possess a hypothetical IR 
attractive “strong-coupling” hxed point, not “visible” within any kind of perturbation 
theory. Then, for d = 1, the huctuation-dissipation theorem along with the Galilean 
symmetry gives the exact values x = 1/2, 2 ; = 3/2 [121 ES]. With additional (rather 
nontrivial) assumptions, one derives dehnite exact values for the exponents in d = 2 
and d = 3 m- Evidence of the existence of the strong-coupling point, provided by 
the so-called functional (also referred to as “exact” or “nonperturbative”) RG [25l |26] , 
although convincing, is numerically still not too impressive, and the situation cannot be 
considered satisfactory; some other open problems are discussed e.g. in [27] 128]. 

It is well known that the behaviour of real systems near their critical points is 
very sensitive to external disturbances, gravity, hnite-size ehects, presence of impurities 
and so on; see e.g. [29l |30] for general discussion and references. What is more, some 
disturbances (randomly distributed impurities in magnets and turbulent mixing of huid 
systems) can change the type of the phase transition or give rise to new universality 
classes with rich and rather exotic properties. 

Investigation of the effects of various kinds of deterministic or chaotic flows (laminar 
shear flows, turbulent convection and so on) on the behaviour of the critical systems (like 
liquid crystals or binary mixtures near their consolation points) has shown that the flow 
can destroy the usual critical behaviour: it can change to the mean-field behaviour or. 
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under some conditions, to a more complex behaviour described by new non-equilibrium 
universality classes [3T] [37]. 

In this paper we study the influence of the random (turbulent) motion of the fluid, 
containing dissolved particles, on the IR behaviour of the randomly growing interface, 
paying special attention to the effects of compressibility. The advection by the velocity 
held v{x) = {nj(x)} is introduced by the “minimal” replacement 

dth Vth = dth + {vidi)h, (1.4) 

where Vt is the Galilean covariant (Lagrangean) derivative. 

We are going to acquire preliminary qualitative understanding of what can happen 
if the fluid motion is taken into account. For this reason, we neglect possible influence of 
the held h{x) on the dynamics of the fluid (“passive” advection) and model the velocity 
held by simple Gaussian statistics with zero mean and prescribed pair covariance with 
vanishing correlation time: 

{Vi{t, x.)vj{t', x')) = 6{t- t') Aj(x - x'), 

/* t/k 1 

Dijiv) = Bo {^p(k) + aQij{k)} exp(ik ■ r), (1.5) 

known as the Kazantzev-Kraichnan ensemble; see e.g. [38] • Here T’jj(k) = 6ij — kikj/k'^ 
and Qij(k) = kikj/k"^ are the transverse and the longitudinal projectors, respectively, 
fc = |k| is the wave number, Rq > 0 is an amplitude factor and a > 0 is an arbitrary 
parameter. The case a = 0 corresponds to the incompressible fluid {diVi = 0), while 
the limit a —)■ cxd at fixed aBo corresponds to the purely potential velocity held. The 
exponent 0 < .^ < 2 is a free parameter which can be viewed as a kind of Holder 
exponent, which measures “roughness” of the velocity held; the “Kolmogorov” value is 
f = 4/3, while the “Batchelor” limit ^ 2 corresponds to smooth velocity. The cutoff 

in the integral (II. bh from below at k = m, where m = l/£ is the reciprocal of the 
integral turbulence scale £, provides IR regularization. Its precise form is unimportant; 
the sharp cutoff is the simplest choice for the practical calculations. 

This ensemble, although it looks simple, has attracted enormous attention in 
turbulence studies because of the deep insight it offers into the origin of intermittency 
and anomalous multiscaling in turbulent advection and turbulence on the whole; see the 
review paper [38] and references therein. The RG approach to that problem is reviewed 
in [39]. In the context of our study, it is especially important that the Kazantzev- 
Kraichnan ensemble allows one to easily model compressibility, which appears rather 
difficult if the velocity field is described by Navier-Stokes equations; see e.g. [IQI E]. 
For a compressible fluid {diVi ^ 0), the covariant derivative can also be introduced in 
an alternative way, namely, = dth -|- di{vih), which is obligatory if the held h has 
the meaning of the density of some conserved quantity. In our case, however, h is “not 
conserved” due to the nonlinear term in fll.2p , and in the following we will consider only 
the variant (II.4p because it preserves the symmetry h h-|-const of the original KPZ 
problem. 
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The plan of the paper is the following. In section [2] we present the held theoretic 
formnlation of the fnll stochastic problem fll.2p . fll.Sp . fll.Sp and diagrammatic techniqne. 
In section [3] we analyze nltraviolet (UV) divergences of the model and demonstrate its 
mnltiplicative renormalizability. Then the RG eqnations, as well as eqnations of critical 
scaling, can be derived in a standard way (section HP . The practical calcnlation of the 
renormalization constants and the RG fnnctions is discnssed in the Appendix. Fixed 
points of the RG eqnations and possible scaling regimes are stndied in section [S] It tnrns 
ont that, in addition to Ganssian hxed point (free held theory), pnrely “kinematic” 
regime (the KPZ nonlinearity is irrelevant in the sense of Wilson) and pnrely KPZ hxed 
point (tnrbnlent transfer is irrelevant), the RG eqnations possess a fnlly nontrivial hxed 
point, in which both the nonlinearity and the mixing are important. The corresponding 
critical exponents can be calcnlated in the form of donble expansions in f and e = 2 —d] 
they are derived in the leading one-loop order. 

The regions of IR stability of the hxed points in the parameter space £, f and a are 
fonnd. In particnlar it tnrns ont, that for small a and most realistic valnes d = 1 or 2 
and .^ = 4/3 or 2, the IR asymptotic behavionr is governed by the kinematic hxed point 
with exactly known exponents. As the degree of compressibility a grows, the stability 
region of the fnll-scale point is getting wider and hnally absorbs the realistic valnes of e 
and f. 

Derivation of the expressions like (II.Ih reqnires (rather simple) analysis of the 
composite helds h^{x)] this is discnssed in section [6l 

It should be admitted, however, that all of our practical results are derived within 
the framework of a standard “perturbative” held-theoretic RG. The strong-coupling IR 
attractive RG hxed point, if it indeed exists, dehnitely survives in the full-scale model, 
but the issue of its IR stability lies far beyound the scope of our study. This problem, 
along with some others, is discussed in sec. [3 


2. Field theoretic formulation of the model 


Let us consider hrst the original KPZ model without the advection. According to the 
general statement [12] (see also the monographs [ini[n]), the stochastic problem (II.2p . 
fll.3p is equivalent to the held theoretic model of the doubled set of helds <F = {h, h'} 
with the action functional 


5 ($) 


-h'Doh' + h' 
2 “ 


—dth + + -(dh)^ 


( 2 . 1 ) 


(we have set Aq = 1). Here and below, all needed integrations over x = (t, x) and 
summations over repeated tensor indices are implied, e.g.. 


h'Doh' = Do 


dt 


dx. h'{t, x) h'{t, x). 


( 2 . 2 ) 


The held theoretic formulation means that various correlation functions and 
response functions of the stochastic problem fll.2p . fll.3p can be identihed with various 
Green functions of the held theoretic model with the action fl2.ip . In other words, they 
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are represented by functional averages over the full set of fields = {h, h'} with the 
weight expiS(<h). 

The bare propagators in the corresponding Feynman diagrammatic techniques are 
determined by the free (bilinear in the helds) part of the action fl2.ip . In the frequency- 
momentum (cj-k) representation they have the forms: 


{hh')o = {h'h)l = 
Do 


{hh)o = 


-io; + ’ 

{h'h')o = 0. 


(2,3) 


+ x^/c^’ 

The model has only one interaction vertex h'{dhY/2. 

In the diagrammatic representation, we will denote {hh)o as a straight line and 
{hh')o as a straight line with a small stroke that corresponds to the held h'. 

Coupling with the velocity held v{x) = {^*( 0 :)} is introduced by the substitution 
fll.dp in fll.21) and thus in fl2.ip . The full problem is then equivalent to the held theoretic 
model of the three helds $ = {h, h' , n} with the action functional 


5 ($) 


]^h'Doh' + h' 


V(/i + + / I" + ^v 


(2.4) 


The last term corresponds to the Gaussian averaging over the field v with the correlator 

(IL3: 

S-^ =-^ J dt J dx j dx'vi{t,x)D~j^{yi-x')vj{t,x'), (2.5) 

where DD is the inverse to the integral operation Dij from fll.5p . 

Thus the Feynman diagrams for the full model fl2.4p involve, in addition to fl2.3p . 
the new propagator fll.5l) and the new vertex —h'{vd)h. 

The role of the coupling constants in the ordinary perturbation theory is played by 
the two parameters 


do ^ V, Wo ^ So/sco ^ 


( 2 . 6 ) 


The last relations follow from the dimensionality considerations (see the next section) 
and dehne the typical UV momentum scale A. 


3. UV divergences and renormalization 

It is well known that the analysis of UV divergences is based on the analysis of the 
canonical dimensions (“power counting”); see, e.g., mm- The dynamic models of 
the type fl2.4p have two independent scales: the time scale T and the length scale L (as 
opposed to single-scale static models). 

Thus the canonical dimension of some quantity F (a held or a parameter in the 
action functional) can be completely descibed by two numbers, the frequency dimension 
df and the momentum dimension d^: 

[F] ~ 
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Table 1. Canonical dimensions of the fields and parameters in the model p.4ll. 


F 

h 

h' 

V 

Xo, X 

^0 

9o 

Bo 

Wo 

g,w,a 

m, /i, A 

dp 

1 

-1 

1 

1 

3 

0 

1 

0 

0 

0 

dp 

-2 

d + 2 

-1 

-2 

-d-4 

2 — d = e 

-2 + e 

e 

0 

1 

dp 

0 

d 

1 

0 

e 

e 



0 

1 


They are found from the obvious normalization conditions 

4 = = 1, d- = d- = 0, dt = = t),dZ = -df = 1, 

and from the requirement that each term of the action functional be dimensionless (with 
respect to the momentum and frequency dimensions separately). Then, based on dp 
and dfp, one can introduce the total canonical dimension dp = dp + 2d‘f (in the free 
theory, dt oc d"^). In the theory of renormalization of dynamical models this factor plays 
the same role as the conventional (momentum) dimension does in static problems; see 
Chap. 5 of [TT] . 

Canonical dimensions of the helds and parameters in the model (12.dh are presented 
in table [H It also includes renormalized parameters (the ones without the subcript “o”) 
and the renormalization mass /i that will be inroduced later on. 

From table [T] it follows the model is logarithmic at d = 2 and ^ = 0, when the 
both coupling constants go and wo simultaneously become dimensionless. Hence, the 
UV divergences in the Green functions manifest themselves as poles in e = 2 — d, f and, 
in general, in all their linear combinations. 

The total canonical dimension of an arbitrary 1-irreducible Green function T = 
($ ... with d) = {h, h', n} in the frequency-momentum representation is given by 

the relation: 

dr = d + 2 — dhNh — dh'Nh' — d^Ny, (3.1) 

where N^, Ny, Ny are the numbers of corresponding helds entering into the function T; 
see, e.g., [H]. 

The total dimension dp in the logarithmic theory (i.e., at £ = ^ = 0) is the formal 
index of the UV divergence: dr = dr|£=.c=o- The superhcial UV divergences, whose 
removal requires counterterms, can be present only in those functions T for which dp is 
a non-negative integer. The counterterm is a polynomial in frequencies and momenta 
of degree dp (provided the convention that a; oc is implied). 

If, for some reason, a number of external momenta occurs as an overall factor in all 
diagrams of a certain Green function, the real index of divergence dp will be smaller than 
dp by the corresponding number of unities. This is exactly what happens in our model: 
the held h enters the vertices h'{dhY and h'{yd)h only in the form of spatial derivatives. 
Thus any appearance of h in some function T gives such an external momentum, and 
the real index of divergence is given by the expression dp = dp — N^. Furthermore, h 
can appear in the corresponding counterterm only in the form of derivative. 




















Random interface growth in random environment 


From table [T] and the expression fl3.ll) one obtains: 

S'^ = 6r-Nh = 4-N,,-- N,. (3.2) 

In dynamical models, all the 1-irredncible Green functions without the response 
helds vanish identically (their diagrams always involve closed circuits of retarded lines); 
see, e.g., m- The sample diagram is shown in Fig. [1] Thus in fl3.2p it is sufficient to 
consider the case Nh/ > 0. 


Figure 1. A diagram of the vanishing function {hh)i-ir with a closed circuit of two 
retarded propagators. 

Then straightforward analysis of the expression (13. 2 h shows that superficial UV 
divergences can be present only in the following 1-irreducible functions: 

{h'h')i_ir ((5r = 0,5p = 0) with the counterterm h'h', 

{h'hh)i_ir ((5r = 2, 5p = 0) with the counterterm h'{dhY, 

{h'h)i_ir ((5r = 2,dp = l) with the counterterm h'd^h, 

{h'hv)i_ir (5r = l,5r = 0) with the counterterm h'{vd)h, 

{h')i-ir (^r = 2, 5p = 2) with the counterterm h', 

{h'v)i-ir (5r = l,5r = l) with the counterterm h\dv), 

{h'vv)i-ir ((5r = 0, (5p = 0) with the counterterm h'v‘^. (3.3) 

Some additional considerations reduce the number of the counterterms. 

The action of the KPZ model is invariant with respect to the transformation 

h{t, x) h{t, X + ut) — It ■ X, h'{t, x) —)■ h'{t, x + ut) (3.4) 

with an arbitrary constant parameter u. This invariance, which becomes the Galilean 
symmetry in terms of the vector field dih, is violated in the full model fl2.ip . However, 
the latter possesses another kind of the Galilean symmetry, namely, 

h(t, x) —)■ h(t, X + ut)^ h\t, x) —)■ h'(t, x + ut), 

n(f, x) —)■ n(f, X + lit) — It (3.5) 

(it is important here that our velocity field is not correlated in time). This symmetry 
requires that the monomial h'{yd)h enter the counterterms only in the form of invariant 
combination h'Vth = h'dth + h'{vd)h. The first term, however, is forbidden by the real 
index (13.2p : the field h appears in it without the spatial derivative. Thus the second 
term is also forbidden (the cancellation of divergent terms from different diagrams can 
be checked in practical calculation, too). The Galilean symmetry also rules out the 
monomial h'v"^. 

The counterterm oc h', stemming from the function {h')i_ir, in terms of the original 
stochastic problem renormalizes the mean value of the random noise (/). This is 
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illustrated by the “tadpole” diagram in Fig. [21 it represents the one-loop contribution to 
the mean value {dhdh). As already mentioned (see the footnote on p. 2), these two mean 
values can simultaneously be ignored. In this respect, the contribution from {h')i-ir is 
similar to the shift of critical temperature in models of equilibrium critical behaviour. 



Figure 2. The one-loop “tadpole” diagram from the function (/i')i-zr- 


The counterterm of the form h'diVi, stemming from the Green function {h'v)i-ir, 
also requires special discussion. It vanishes identically for the incompressible case, where 
diVi = 0. However, the practical one-loop calculation shows that it is absent in the 
general case {a ^ 0). One can give some arguments that this is true in all orders of 
perturbation theory; see section Appendix A Since our present analysis is restricted 
to the one-loop calculations, we will not take this term into account in the following 
discussion. 

Then we are left with the three counterterms of the form h'h', h'd^h and h'{dhy. 
All these terms are present in the action (12.dh making our model multiplicatively 
renormalizable. The renormalized action then can be written in the form: 

5r($) = \zih'Dh' + h' {-dth + KZ2d‘^h + \z:,{dhf + / I + (3.6) 


Here g, w and x are renormalized analogs of the bare parameters, and the functional 
Sv from fl2.5p should also be expressed in renormalized variables. The renormalization 
constants Zj depend only on the completely dimensionless parameters g,w,>i,a and 
absorb the poles in £ and f. 

The renormalized action (13.6p is obtained from the original one (12.dh by the 
renormalization of the helds h —)■ Zhh and h' —?• Zh'h' and of the parametrs: 


Kq 


g^ gja Zg^ wq w^^Z^. 


(3.7) 


The amplitudes D and B are expressed in renormalized parameters as follows: 
D = gx^gf, B = 


(3.8) 


The renormalization constants in Eqs. fl3.6p and fl3.7p are related as follows: 

Zg = Z1Z2 Z^ = Z2, Zh = Z3, Z^i = Z3 

ZhZh' = 15 = 1, ZyjZ^ = 1. (3.9) 


The hrst two relations in the second line follow from the absense of the counterterm 
h'Vfh, the last relation follows from the absense of renormalization of the term S^- 

The renormalization constants Z 1 -Z 3 are calculated directly from the diagrams, 
then the constants in fl3.7p are found from fl3.9p . The renormalization constants can be 
found from the requirement that the Green functions of the renormalized model (13.61) 
be UV hnite when expressed in renormalized variables. In our case this means that the 
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Green functions are UV finite at ^ ^ 0. The calculation in the hrst order in g and w 
(one-loop approximation) gives (see section Appendix A for the details): 


Zi = 1- q - -w, 


Zo — Z:t — 1 — 


w d — 1 + a 




2d 


(3.10) 


where g = gSd/{27iY, w = wSd/{27iY, and Sd = 27i‘^/r{d/2) is the area of the unit 
sphere in d dimensions. 

If the minimal subtraction (MS) scheme is employed, the renormalization constants 
must have the forms “Z = 1 -|- only poles in e and (and in higher orders in their 
linear combinations). Then, strictly speaking, in our one-loop accuracy we have to 
replace d = 2 — e ^2 in the above expressions. However, for some time we will 
keep them in the form (13.10^ : then some exact results for the special cases will be 
derived; see section [5l A similar renormalization scheme, where the dimension d is kept 
in “geometrical factors” stemming from contractions of various projectors, was earlier 
used in [ 21 ]; its validity and equivalence to the MS scheme was demonstrated in [ 22 ] . 


4. RG equations and RG functions 

Consider briefly an elementary derivation of the RG equations; detailed exposition can 
be found in the monographs [TillITT] . The RG equations are written for the renormalized 
Green functions Gr = (<!)■■■ <h)R. They differ from the original (unrenormalized) 
Green functions G by overall numerical factors (due to rescaling of the helds) and 
by different choice of the parameters (e,/i instead of cq). Thus the renormalized 
Green functions can be equally used for analyzing the critical behaviour. The relation 
^^(Z^d), e,/i) = S'(<h,eo) between the functionals ( 12 .dh and (13.61) yields the relations 

G(e„,,,,) = Z""z"«'GR(e,^,,,,), (4,1) 

between the Green functions. Here, as above, Nh and N^i are the numbers of 
corresponding helds entering into G (we take into account that in our model Z^ = 1 ); 
Co = {do, ^ 0 , is a full set of bare parameters and e = {g, x, w} are their renormalized 
counterparts; the ellipsis stands for the other arguments (times, coordinates, momenta 
etc.). 

We use to denote the differential operation /r9^|eo- When expressed in the 
renormalized variables it looks as folows: 

T^rg = + l^gdg + jdwdw — (4-2) 

where Vx = xdx for any variable x. The anomalous dimensions 7 are dehned as 

7 ^ = V^inZp for any quantity F, (4.3) 

and the ft functions for the two dimensionless couplings g and w are 

fdg = V^g = g[-e- 7 g], (3^ = V^w = w[-f- 7 ^], 

where the second equalities come from the dehnitions and the relations (13.7p . 


(4.4) 
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In order to derive the basic RG differential equations we apply to the both sides of 
the equality fl4.ip the operation V^: 

{T^rg + Nhlh + Nh'^lh'} GR{e, /r,...) = 0. (4.5) 

At last, equations fl3.9p yield the following relations between the anomalous 
dimensions fl4.3p : 

7h = 73, 7/*' = -73, Iw = -7x, 'yv = 1, 

7x = 72, 7g = 7i - 372 + 273 . (4.6) 

The anomalous dimension corresponding to a given renormalization constant Zp 
can be found from the relation 


7 f = {l3gdg + jSnjd^) In Zf ~ - {eVg + In Zp, 


(4.7) 


obtained from the definition fl4.3p . expression fl4.2l) for the operation in renormalized 
variables, and the fact that the renormalization constants depend only on the two 
completely dimensionless coupling constants g and w. In the second part of the relation, 
we retained only the leading-order terms in the (3 functions fl4.4p as it is sufficient for 
the first-order approximation. Utilizing the MS scheme, in the one-loop approximation 
from the explicit expressions fl3.10p one finds: 


'yi = g/% +aw/2, 72 = 73 = u) 


d — 1. a 
2d ’ 


(4,8) 


with g and w defined earlier (I3.10p and the corrections of order g‘^, wg and higher. 


5. Fixed points and scaling regimes 

It is well known that the long-time large-distance asymptotic behaviour of a 
renormalizable field theory is determined by IR attractive fixed points of the 
corresponding RG equations. In general, the coordinates of possible fixed points are 
found from the requirement that all the (3 functions vanish. In the model fl2.4p the 
coordinates g^, are determined by the two equations 

/3g{g*,w^) = 0, /3n,{g*,w^) = 0, (5.1) 

with the (3 functions given in fl4.4p . The type of a fixed point is determined by the 
matrix 

G = = d^i/dgj}, (5.2) 

where /3i is the full set of the (3 functions and gj = {g,w} is the full set of coupling 
constants. For an IR attractive fixed point the matrix G should be positive, i.e., the 
real parts of all its eigenvalues should be positive. 

From the relations (14.6p we find that 7g = 71 — 72, 7«, = —7x = 72, so that 


1^9 = g (-^ -lg)= g (-^ - 7i + 72), 

I 3 w = w (-f - 7^) = w (-^ - 72). 


(5.3) 
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Substituting explicit expressions fl4.8p we arrive at the explicit one-loop expressions for 
the ft functions: 


Pg = 9 
ftw = W 


q ^ id — l){a 

8+”^- k- 


-f + w 


{d — 1 -|- cr) 
2h 



(5.4) 


with possible corrections of the order w'^ and so on. It is worth noting that the one-loop 
result fl5.4p for (3^ becomes exact at = 0 (see e.g. [31]), while the result for ftg becomes 
exact at tc = 0 (as follows from the analysis of |23]). 

From Eqs. fl5.ip and fl5.3p one hnds that there are four different hxed points in our 
model. As dgftw = 0, the matrix is triangular in every case and its eigenvalues are 
simply given by the diagonal elements fig = dftg/dg and fl^ = dftw/dw. 

The hxed points are as follows: 


1. Gaussian (free) hxed point: g^, = = 0; fig = —flw = —f (all these expressions 

are exact). 

2. w* = 0 (exact result to all orders), g^ = —Se; fig = £, fl^ = —f- 


This hxed points corresponds to pure KPZ model: although the interaction with 
the velocity held is present, it does not ahect the leading-order IR asymptotic behaviour 
(it is irrelevant in the sense of Wilson). 


3. (?* = 0 (exact), = 


2df 


(d — 1 -|- cr) 


] fig — —e + f — 


da f 


{d — 1 -|- o) 


, f (exact). 


This hxed points corresponds to pure Kraichnan model with small-scale stirring; 
the KPZ nonlinearity is irrelevant in the sense of Wilson. 


4. g* = 

(exact). 


-e + i- 


daf 

{d — 1 -|- of 


w. = 


2df 


(d — 1 -|- cr) 


] fig — e — f-\- 


da f 


(d — 1-1- o) 


, f^w ^ 


This hxed point corresponds to a new nontrivial IR scaling regime (universality 
class), in which the nonlinearity of the model fl2.4p and the turbulent mixing are 
simultaneously important. 

In hgure [3] the regions of IR stability for all the hxed points in the e-f plane are 
shown: dark space for the Gaussian point, horizontal shading for the KPZ point, vertical 
shading for the Kraichnan point and white space for the new regime. These regions are 
areas in which the eigenvalues of the matrix fl5.2p for the given hxed point are both 
positive. 

In the one-loop approximation (15.3p . all the boundaries of the regions of stability 
are straight rays. Diherent regions have neither gaps nor overlaps between them. Such 















Random interface growth in random environment 


13 


A 




Figure 3. Regions of stability of the fixed points in the model (I2.4I1 . 


a pattern is typical feature of the first-order approximations. The boundaries £ < 0, 
.^ < 0 for point 1, £ > 0 for point 2 and > 0 for point 3 are exact, while the others 
can be affected by the higher-order corrections, i.e., the boundaries may become curved 
and gaps or overlaps may appear between the different regions of IR stability. 

The main qualitative conclusion that can be drawn from this pattern, is that for 
small a and most realistic values d = 1 or 2 and ^ = 4/3 or 2, the IR asymptotic 
behaviour is governed by the Kraichnan hxed point. However, as the degree of 
compressibility a increases, the stability region of the new point is getting wider and 
hnally absorbs the realistic values of £ and f. Indeed, the boundary between the regions 
3 and 4 depend on a. When a grows, it rotates counterclockwise and, for a —)■ cx), 
approaches the ray e = {1 — d)f. Thus for d = 1 it tends to the vertical ray £ = 0, ^ > 0 
and for d = 2 it tends to £ = —f, ^ > 0; see hgurelH (note, that for d = 2 the boundary 
becomes vertical at a = 1). 

We recall, however, that the results concerning the boundary between the regions 
3 and 4 may be affected by the higher-order contributions. We also note that 
interpretation of the hxed point 4 leads to the same difficulty as that of the KPZ point: 
when it is IR attractive, the coordinate lies in the unphysical region < 0. We will 
return to this issue in section [71 

6. Critical dimensions. Critical scaling of the structure functions. 
Composite fields h^{x). 

Existence of an IR attractive hxed point implies existence of scaling (self-similar) 
behaviour of the Green functions in the IR range. In this critical scaling all the 
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/ 
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Figure 4. The boundary between the regions 3 and 4 moves counterclockwise with 
the growth of the compressibility a until it reaches the ray e = —^ (for d = 2). 


“IR irrelevant” parameters (/r, g and w in onr case) are kept fixed, while the “IR 
relevant” parameters (coordinates/momenta, times/freqnencies, the helds) are dilated. 
The critical dimension Ap of a certain IR relevant qnantity F is given by the relations 
(with the normalization condition = 1 ) 

Ap = dp + Ai^dfp + 'jp, ( 6 . 1 ) 

where 

A-= 2-7: (6.2) 

is the critical dimension of freqnency, d’f.‘^ are the canonical dimensions of F, given in 
table [H and 7 ^ is the valne of the anomalons dimension (14.311 at the fixed point in 
qnestion: 7 ^ = see e.g. [H] for detailed explanation. 

From table [1] for the dimensions of the fields we obtain 

Ah = —2 + A^ + 7 ^, Ah' = (d + 2) — + 7 )). 

Then the relations (14. 6 h give 

A^ = 2 - 7 *, A;, = - 7 ; + 7 *, Ah' = d + 7 ; - 7 *. 

Finally, snbstitnting explicit one-loop expressions (14. 8 h yields 

Ah = 0, Ah' = d, Aoj = 2 (6.3) 

for the fixed points 1 and 2 and 

A;, = 0, Ah' = d, A^ = 2-f (6.4) 

for the points 3 and 4. For the points 1-3 these expressions are exact; see the remark 
below eq. (15.411 . For the point 4 the dimensions of the fields can be affected by the higher- 

order contribntions with the exact condition that Ah + Ah' = d, and the expression for 

A^ is exact dne to the second relation in (15.311 . 

In the traditional notation (II.Ih one has Ah = —x and A^^ = 2 ;. However, the 
qnantity Sn in (11.111 is not an ordinary n-th order Green fnnction of the basic helds h{x): 
it is a snm of pair correlation fnnctions {h^{x)h'^{0)) of the composite helds (“composite 
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operators” in the quantum-field terminology) h'^{x). Renormalization of such quantities 
requires further analysis which, however, is rather simple in the present case. 

Total canonical dimension of the 1-irreducible Green function T = (F<h... 
with one composite operator F and arbitrary number of basic fields <h = {h^h',v} in 
our model is dr = dp — dhNh — dh'Nh> — d^Ny, where Nh, Nh', Ny are the numbers of 
corresponding fields entering into the function T, dh,h',v are their canonical dimensions 
and dp is the canonical dimension of F; see, e.g., m- The formal index of divergence is 
dr = dr|£=^=o; superficial divergences can be present in T if dr is a non-negative integer. 

From table [1] for F = h" we obtain dp = 0 and dr = — Nv- Thus the 

divergences, at first sight, can be present in all functions T = ifi^ h... h)i-ir with 
Nh' = Ny = t], arbitrary and dr = 0. However, any nontrivial diagram of such 
a function involves at least one external vertex h'{dhY or h'{vd)h, where the field h 
stands under a derivative. Thus at least external momentum appears in the diagram as 
an overall factor, the real index dp is negative and the superficial divergence is in fact 
absent. 

This means that all the operators F = are not renormalized and their critical 
dimensions are simply given by AiT’ = nAh- This justifies the relation (II.ip with the 
dimensions fl6.3p . fl6.4p . 

7. Discussion and conclnsion 

We studied effects of turbulent mixing in the problem of randomly growing interface. 
The growth was modelled by the well-known Kardar-Parisi-Zhang stochastic equation 
fll.2p . fll.Sp . The turbulent velocity field was modelled by the Kraichnan’s rapid-change 
ensemble fll.Sp . 

The full problem can be refornulated as the multiplicatively renormalizable model 
with the action functional fl2.4p . Using the field theoretic RG we show that, depending 
on the relation between the exponent f and the spatial dimension d, the system exhibit 
different types of IR behaviour, associated with four possible fixed points of the RG 
equations. In addition to known regimes (ordinary diffusion, ordinary growth process, 
and passively advected scalar field), existence of a new nonequilibrium universality class 
is established. 

Practical calculations of the fixed point coordinates, their regions of stability and 
critical dimensions were performed in the first order of the double expansion in f and 
e = 2 — d (one-loop approximation). 

It was shown that for the incompressible fluid the most realistic values of f and d 
correspond to the universality class of passive scalar field, when the nonlinearity of the 
KPZ model is irrelevant. If the degree of compressibility a becomes large enough, the 
crossover in the IR behaviour occurs, and these values of d and f fall into the region of 
stability of the new regime. 

However, a few problems remain open. 

For the new universality class (as well as for the KPZ model) the coordinates 
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of the fixed point lie in the unphysical region 9* < 0 , which corresponds to the 
“wrong” negative sign of the amplitude in the pair correlator fll.Sp . Thus it requires a 
careful physical interpretation. In this connection one can recall, that in the Doi-Peliti 
formalism jHJ |l5] , where the original microscopic problem is formulated in terms of the 
creation-annihilation operators, the terms quadratic in the response helds can appear 
in the action functionals with the negative signs; see e.g. |45] . 

Another question is the fate of the strong-coupling hxed point of the pure KPZ 
model [25]. If it indeed exists, it dehnitely survives in our problem (with the second 
coordinate w* = 0), but can become unstable. On the other hand, new nonperturbative 
hxed points with tc* 7 ^ 0 can appear. 

In our analysis we employed the simplest Kraichnan’s ensemble for the advecting 
velocity held. It would be interesting to consider more realistic models with non- 
Gaussianity, hnite correlation time, anisotropy, etc. This work is in progress. 
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Appendix A. Calculation of the renormalization constants 

In this Appendix we shortly present derivation of the hrst-order results (lO.lOh for the 
renormalization constants. Although the one-loop calculation is rather simple and can 
be accomplished in a few diherent ways, it is worth to discuss it for completeness and 
in order to mention some interesting subtleties specihc of the model fl2.4p . 

The renormalization constants are determined by the requirement that the Green 
functions of the renormalized model (13.bh . expressed in renormalized variables, be UV 
hnite (in our case, be hnite at £ —)■ 0, 0). The full set of constants Zi-Z^ can be 

found from the three 1 -irreducible functions: {h'h)i-ir, (h'h')i_ir and {h'hh)i-[r- In the 
renormalized model, the corresponding one-loop approximations have the forms 



(A.l) 



(A.2) 



and 
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(A.3) 


Here we denoted the bare propagator (hh)o as a straight line, (hh')o as a straight 
line with a small stroke that corresponds to the field h', and the velocity propagator as 
the wavy line. In the following, we set the external frequency rj equal to 0, because the 
divergent part of the function (h'h)i_ir is proportional to where p is the external 
momentum. 

All the diagrammatic elements should be expressed in renormalized variables using 
the relations fl3.6p - fl3.9p . In the one-loop approximation, the constants Zi in the bare 
terms of flA.ip . flA.2p and (IA.3p should be taken in the first order in g and w, while in 
the one-loop contributions they should simply be replaced with unities, Zj —)■ 1. Thus 
the passage to renormalized variables in the one-loop diagrams is achieved by the simple 
substitutions xq —>■ x, po ^ and wq wpf. 

The IR regularization in the diagrams involving the velocity propagator is provided 
by the cutoff in the integral (11.5p from below at k = m. In other diagrams, IR 
regularization is provided by external momenta and frequencies. We are interested, 
however, only in the UV divergent parts of these diagrams (poles in e and f). Thus 
we will use the following trick, which simplifies the calculation: integrations over the 
momenta in all diagrams will be cut off from below at k = m. Then in the logarithmically 
divergent functions flA.2p and flA.311 external momenta and frequencies can be set equal 
to zero, when in the quadratically divergent function flA.ip we will keep only the p^ term 
of the expansion in the external momentum p. 

Now the integrations over the frequency are easily performed by residues. The 
resulting integrals over the momentum with the aid of the formulas 

jdkkif{k) = 0, jdk^^f{k) = ^ / 


where f{k) is any function depending only on k = |k|, are reduced to the scalar integral 

1 


J{m) = 


dk 


' k>m 


kd+y 


= S, 


m 


-y 


y 


(A.5) 


with Sd from (I3.10p . Here either y = e or y = f. 

The last two diagrams in flA.3p in fact vanish and thus give no contribution to the 
renormalization constant. Indeed, they effectively involve closed circuits of retarded 
propagators; it is crucial here that the velocity correlator contains the 6 function in 
time. 

Direct calculation shows that the first three diagrams in (1A.3p also give no 
contribution to because their divergent parts cancel each other. This is a consequence 
of the Galilean symmetry (13.4p of the original KPZ model, which forbids the counterterm 
h'{dhy in all orders of perturbation theory. 
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The analytic expression for the only remaining diagram in flA.Sp has the form: 


du 


dk 




J W) (W w 

Here the prefactor comes from the two lower vertices, the hrst cofactor in the 
integrand comes from the npper vertex (nnmerator) and from the propagators {h'h)^ 
(denominator), and the remaining factor is the velocity correlator. Proceeding as 
explained above, we hnally obtain for flA.3p : 

/r \ ^ {d — 1 + a) 


{h'hh)i_ij. — < Z3 + — ^ 


f \m/ 


2d 


(A.7) 


The factor is UV hnite: it tends to nnity for 0. We can see that, in order 

to cancel the pole in f in (IA.7p . the renormalization constant Z 3 can indeed be chosen 
in the form (I3.10p . 

Now let ns tnrn to the renormalization constant Z 2 . 

The hrst one of the two diagrams in flA.lD appears UV hnite and does not contribnte 
to Z 2 . To see this, consider the corresponding analytic expression (np to insnfhcient 
amplitnde factors): 


duj 

W) 


dk ki{p + k)ikjPj 


lk>m (2vr)^ 

dk ki{p + k)ikjPj 


-iuj + x|k + p| 


oc 


y,>^(2vr)'^F(F+|k + p|2)- ^ ^ 

We are interested in the p^ term of its expansion in p. It is a snm of two identical 
integrals with opposite signs. Indeed, the hrst one comes from the contribntion piPjkikj 
in the nnmerator; then in the denominator we can pnt p = 0. The resnlting integrand 
becomes eqnal to PiPjkikj/2k‘^. The second one comes from the contribntion kikikjPj. 
Then one has to expand the denominator np to the order 0(p): 


dk 




dk kjPj 


lk>m (2vr)'^ P + |k + p| 


' k>m 


{2'kY 2k‘^ 


1 - 


(Pk)\ 

j 


(A.9) 


The hrst term in flA.9D vanishes becanse the integrand is odd in k, and the second eqnals 
to the aforementioned one np to the minns sign. 

The analytic expression for the second diagram in flA.ip is as follows: 


du 

W) 


' k>m 


dk ipiiip + k)j . {-Pp(k) + aQij{k)} 

Wf ~k^ + ■ 


Integration over u involves the indeterminacy 

^ ' = 9(0). 


(A. 10) 


(A. 11 ) 


J (27r) —ia; + x|k + p| 

where d(0) is the Heaviside step fnnction at the origin. This rehects the details of the 
velocity statistics lost in the white-noise limit; see the discnssion in |38]. In the case 
at hand, the 6 fnnction in fll.Sp shonld be nnderstood as the limit of a narrow fnnction 
which is necessarily symmetric in t, T, becanse one deals with a pair correlator. Thns 
the indeterminacy in flA.llD must be unambiguously resolved as half the sum of the 
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limits: 6*(0) = 1/2. The integrand of the resulting integral over k in (lA.lOD has an odd 
contribution that can be dropped and the even contribution that yields 


(h'h)i_ir = 


V , w f {d-l +a) 

+ “ — Td — 


(A.12) 


f \mJ 

We can see that Z 2 can be chosen in the form fld.lOp : in the one-loop approximation 
Z 2 = ^ 3 - 

We will not discuss the calculation of Zi in detail; the first diagram reduces to 
the scalar integral flA.Sp with y = e immediately after integration over the frequency. 
The analytic expression for the second diagram is similar to flA.bp with the replacement 
PiPj —)■ kikj, so that only the contribution from the longitudinal projector survives. This 
leads to the replacement {d — 1 + a)/2d —)■ a. Bringing all contributions together and 
taking into account the symmetry coefficient 1/2 for the first diagram, one obtains: 


8 e \m/ 2t \m/ 


2 ? 


(A.13) 


We can see that the constant Zi that removes the poles from expression (IA.13H can be 
taken in the form fId.lUp . 


Figure Al. The one-loop contribution to the function {h'v)i-ir. 

It remains to discuss possible UV divergence in the function {h'v)i-ir with the 
counterterm h'diVi. The only one-loop diagram for this function is shown in hgnre lATl 
It is not difficult to see that the corresponding analytical expression is nearly identical 
to that of the hrst diagram in the function flA.ip : the latter only has additional extra 
factor pj. Thus the diagram in question is also UV hnite and does not give rise to the 
corresponding counterterm. 
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